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First-principles electronic-structure methods are used to study a structural model for 
Ag/Si(lll)3x 1 recently proposed on the basis of transmission electron diffraction data. The fully 
relaxed geometry for this model is far more energetically favorable than any previously proposed, 
partly due to the unusual formation of a Si double bond in the surface layer. The calculated electronic 
properties of this model are in complete agreement with data from angle- resolved photoemission and 
scanning tunneling microscopy. 
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The surfaces of silicon reconstruct in strikingly diverse 
ways. This diversity provides a rich proving ground 
for simple, physically intuitive ideas about the stabil- 
ity of semiconductor surfaces — ideas which are invalu- 
able for understanding more complex dynamical phenom- 
ena such as growth, etching, and reactivity. Two such 
simple concepts — elimination of surface dangling bonds 
and relief of surface stress — explain the frequent appear- 
ance of elementary "building blocks" in silicon recon- 
structions. For example, dimers appear on Si(001)2xl 
(and its variants), (111)7x7, (113)3x2, (114)2x1, and 
(5,5,12)2x1; 7r-bonded chains appear on (111)2x1 and 
(5,5,12)2x1; and adatoms appear on the (111)7x7, 
(113)3x2, (114)2x1, and (5,5,12)2x1 reconstructions 
(see Ref. [Q and citations therein). 

One reconstruction that remains controversial is the 
metal-induced M /Si(lll)3xl (where M =Li, Na, K, Ag, 
Mg), which is widely believed to have a single common 
structure. Starting from known building blocks, we and 
others have recently proposed two models for this recon- 
struction, both with low dangling-bond density and low 
surface stress |^|-^| . First-principles total-energy calcula- 
tions showed these models to be stable relative to previ- 
ously proposed ones ||. The more stable of the two, the 
extended Pandey chain model, was also consistent with 
scanning tunneling microscopy (STM) images, but the 
calculated surface-state band structure of both models 
was in serious disagreement with angle-resolved photoe- 
mission (ARPES) data §. 

In this Letter, we examine theoretically a new model 
for the M:3 x 1 surface proposed very recently by Collazo- 
Davila, Grozea, and Marks (CGM) on the basis of direct 
phasing of transmission diffraction data || , and indepen- 
dently by Lottermoser et al. from surface x-ray diffraction 
and total-energy calculations First, starting from co- 
ordinates obtained by CGM, we further relax the atomic 
positions so as to minimize the calculated total energy. 
The resulting model is by far the most stable of any pro- 



posed to date. Second, we show how a nearly perfect 
"surface symmetry" of this reconstruction neatly resolves 
a long-standing puzzle regarding ARPES data for Li:3xl 
(and may also explain the apparent insulating nature 
of Mg:3xl, which has an odd number of electrons per 
unit cell) . Third, we show that this model derives its re- 
markable stability from the formation of a true Si double 
bond — a "building block" not seen on any other surface 
of silicon. Finally, we show that this model completely 
accounts for the appearance of the existing STM images 
(including the differences between M=Li and Ag). 

The model proposed by CGM was based on data from 
Ag:3xl J^]. Previous experimental results from low- 
energy electron diffraction Jl(J, STM [0, and core level 
spectroscopy HQ suggest that Li, Na, K, Ag, and Mg all 
have the same reconstruction (apart from several details 
discussed later). For this reason we concentrate first on 
Li:3x 1, which leads to a slightly simpler physical picture, 
and then later discuss the differences between Li:3x 1 and 
Ag:3xl. Starting from the atomic positions determined 
by CGM, we performed full structural relaxation using 
first-principles electronic-structure methods. The calcu- 
lations used a double-sided slab geometry with 14 layers 
of Si plus Li, with a vacuum region equivalent to 10 lay- 
ers of Si. The four innermost layers of Si were fixed, 
while all other atoms were allowed to relax until the rms 
force was less than 0.02 eV/A. Total energies and forces 
were calculated within the local-density approximation 
(LDA), using Troullier-Martins pseudopotentials and a 
plane- wave basis with a kinetic-energy cutoff of 12 Ry, 
as implemented in the fhi96md code [Q . Four fc-points 
were used to sample the surface Brillouin zone. 

The fully optimized structure is shown in Fig. [l]. Sev- 
eral observations will motivate our discussion below. (1) 
There are four incquivalent Si surface atoms (labeled a, 
b, c, d in Fig. |l|) which form a "honeycomb chain" ly- 
ing in a plane parallel to the surface. (2) Each of these 
surface Si atoms is threefold coordinated: outer atoms a 
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FIG. 1. Fully relaxed geometry for the honeycomb 
chain-channel model of Li/Si(lll)3x 1. Black circles are Li. 
White, light gray, and dark gray circles are the surface layer, 
first layer, and deeper Si layers, respectively. The heavy dot- 
ted line is an approximate mirror-symmetry plane for the top 
two Si layers (see discussion in text). The 3x1 surface unit 
cell is shown as a heavy solid line. A disilene molecule, Si2H4, 
is shown in the lower right (see discussion). 

and d roughly tetrahedrally, and inner atoms b and c in 
a planar configuration. (3) Electron density plots show 
that b and c are only very weakly bonded to the first- 
layer Si atom, e, below. (4) The Li atom sits in a channel 
formed by neighboring honeycomb chains, and occupies a 
symmetric location with respect to the [110] direction — 
that is, it preserves the [TlO] mirror plane. (5) Electron 
density plots show that the Li atom is fully ionized. 

To compare the stability of this model with previous 
ones we compute the surface energy, E s — [E t (N) — 
NE b t ulk ]/2, where E t {N) is the total energy of a re- 
constructed slab supercell containing N Si atoms, and 
E\ ulk is the total energy per atom of bulk Si. Relative 
to the extended Pandey model, which is the lowest en- 
ergy model to date ||, this "honeycomb chain-channel" 
(HCC) model has a surface energy lower by 0.52 eV per 
3x1 cell, or 14 meV/A 2 — an enormous energy gain on 
the scale of silicon surface reconstructions (and in good 
agreement with the LDA result of Rcf . Q ) . 

A structural model with five under-coordinated atoms 
would appear to be energetically unfavorable, so how can 
we understand the very low calculated surface energy of 
the HCC model? We develop here a simple intuitive pic- 
ture for this stability, one which also resolves an appar- 
ent discrepancy regarding ARPES data for Li:3xl. The 
measured spectra reveal a single surface state, dispers- 
ing downward from T toward the A point, i.e., along the 
[TlO] direction 0. The apparent absence of additional 
states is puzzling since, roughly speaking, each surface 
orbital is expected to give rise to a surface state — and 
the 3x1 structural models proposed to date have at least 
three surface orbitals. Indeed, the HCC model has five 
unsaturated Si surface orbitals, and so should give rise 
to five surface states. There are six electrons available 
to fill these states — four from the Si surface atoms (a, b, 
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FIG. 2. Theoretical band structure for the honeycomb 
chain-channel model of Li:3xl. Surface states, S, are denoted 
by solid curves: heavy lines for states of even parity, and light 
lines for states of odd parity. The projected 3x1 bulk bands 
are shaded. Experimental angle-resolved photoemission peak 
positions are indicated by circles. 

c, d), one from the Li ion, and one from e — and so we 
expect three filled surface states and two empty states. 

The theoretical band structure, part of which is shown 
in Fig. H confirms this description and illustrates the dis- 
crepancy. Three surface states, labeled S^~, S^~, are 
fully occupied and one, labeled , is completely empty 
(the other empty state lies higher and is not shown). For 
the occupied states, the single peak in the ARPES data, 
also shown in Fig. ^, lies between and thus is consistent 
with and , but is completely inconsistent with . 

The resolution to this discrepancy provides internally 
consistent evidence in support of the HCC model itself, 
and also explains its surprising stability. Our argument 
is that the photoemission cross-section is large for states 
Si and S2 , but small or vanishing for S2 ■ Specifically, 
if the plane defined by the incident light beam and the 
photoelectron detector is a mirror plane of the crystal, 
then the states along that direction in reciprocal space 
have either even or odd parity. Given this condition, the 
photoemission matrix element, i\A ■ p\^! f) , will vanish 
by symmetry for two situations: (i) when the polarization 
vector, A, is normal to the mirror plane and the initial 
wave function, has even parity; or (ii) when A is 
parallel to the mirror plane and has odd parity jl3| . 

The ARPES data shown in Fig. || are from a syn- 
chrotron light source with the incident beam, detector, 
and polarization vector all lying in the plane defined by 
the surface normal and the [TlO] direction, i.e. the plane 
normal to [112]. This plane is shown in Fig. |l| as a heavy 
dotted line. It is clearly not a mirror plane with respect to 
either the Li atoms or the Si substrate. However, the Li 
atoms are completely ionized, and the substrate atoms do 
not contribute significantly to the surface states. Hence 
neither is relevant to determining the surface symmetry. 
The relevant atoms are those in the surface layer (a, b, c, 
d) and perhaps the first layer. Fig. |l| shows that for these 
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FIG. 3. Wave functions (squared) for the four surface states 
shown in Fig. ^, evaluated near k=A. The plotting plane is a 
vertical cut perpendicular to [TlO] and passing through (from 
left to right) surface atoms b, c, a d (projected positions are 
marked below). Large black circles denote atoms in the plot- 
ting plane, and small circles denote atoms out of this plane. 
Contours are logarithmically spaced. 



atoms the dotted line is a nearly perfect mirror plane. 

The extent to which this approximate symmetry leads 
to surface states with definite parity can be easily evalu- 
ated. Given an exact mirror-plane symmetry and a set of 
orbitals, <fi, associated with atoms a through e, the five 
symmetry-adapted basis functions will have even (+) and 
odd (— ) parity: 
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b a ± <t>d: $3 



(1) 



In Fig. [| we show the squared modulus of the four sur- 
face states from Fig. [^, evaluated near the A point along 
the [110] direction. This figure clearly shows that is 
a linear combination of &f and 4>g" and thus has overall 
even parity (the orthogonal state, not shown, is ~2 eV 
above the valence edge), while S± « $^ has odd parity, 
as evidenced by the node at the mirror plane. Likewise, 
S2 ~ $ 2 and S2 &2 have overall even and odd parity, 
respectively. Having shown that the existence of an ap- 
proximate surface mirror-plane symmetry indeed leads to 
surface states with approximate parity, we conclude that 
the matrix element for photoemission from the states S^ 
and S2 will be very small. This explains the absence 
of any signal from S% in the ARPES data. Moreover, 
we note that the single peak in the ARPES spectra is 
quite broad, with a FWHM of 0.7-0.8 eV. We suggest 
that this broad peak is in fact a superposition of states 
and S£- Indeed, excellent fits of the zone boundary 
spectrum can be obtained by using two Voigt functions, 
representing S± and S%, separated by 0.4 eV. 

The description of the occupied surface states in terms 
of a symmetry-adapted basis also leads to a simple expla- 
nation for the stability of the HCC model. The largest 
contribution to the S± state is from the &f basis func- 
tion, which is an even combination of orbitals separated 
by one Si bond length — and thus bonding with respect 



to b and c. The antibonding combination of these or- 
bitals, Si ~ is unoccupied. This suggests the exis- 
tence of a true Si double bond, between atoms b and c, 
which we propose is primarily responsible for the stability 
of the HCC model. Indeed, the energy splitting between 
and is almost 2 eV, far larger than the bonding- 
antibonding splitting (less than 0.5 eV) for the extended 
Pandey chain. 

The Si double bond, although unusual, also occurs in 
a variety of molecular systems. In the first known ex- 
ample, the Si=Si bond was sterically stabilized by two 
mesityl groups bonded to each Si Q . The simpler disi- 
lene molecule, Si2H4, has a planar structure (see Fig. [j]) 
essentially identical to the arrangement of atoms a, b, c, 
d in the HCC reconstruction. Here, the bonds from the 
substrate to a and d play the role of steric stabilizers. 

We turn now to the appearance of the M/Si(lll)3xl 
surface in STM images. Experimental filled-state images 
of Li:3xl appear as double rows of staggered maxima 
with approximately equilateral spacing; empty-state im- 
ages appear as narrow lines with small spurs on both sides 
]rl| . It is at first difficult to reconcile the STM image of 
staggered maxima with the model of Fig. |l|, which has a 
mirror plane at the center of each chain. This, too, has 
a natural explanation. The 3x1 surface periodicity im- 
plies a half-period shift between adjacent chains along the 
[TlO] direction. Thus, orbitals 4> a and <pd from two adja- 
cent chains should give rise to a double row of spots with 
approximately equilateral spacing. Both S2 ~ 4> a + <i>d 
and S2 ~ 4>a — 4>d lie within 0.5 eV of the valence-band 
maximum, so that filled-statc images with this bias volt- 
age or larger should appear as described. 

In Fig. ^ we show theoretically simulated constant- 
current STM images for the HCC model, obtained by 
integrating the local-state density over a 1.5-eV energy 
window for both filled and empty states. The filled-state 
simulated image confirms our qualitative description, and 
is in excellent agreement with experimental STM images 
of Li:3xl. (Note also a small contribution from the S^ 
state, which contributes to this simulated image because 
of the large bias voltage chosen.) The empty-state sim- 
ulated image is also in excellent agreement with experi- 
ment, appearing as a single row, centered on the Li ions, 
with small side spurs. This is easily understood as the 
superposition of the unoccupied Li 2s atomic state and 
the unoccupied S-f surface state. We also remark that 
the registry of the filled- and empty-state images has been 
determined experimentally with respect to both the [TlO] 
and [112] directions Jul, and is in identical to that of the 
simulated images in Fig. [|. 

We return now to the differences between M=Li and 
Ag. We repeated the above procedure, starting from the 
geometry determined by CGM for the Ag:3xl surface 
and using an appropriately higher energy cutoff of 40 Ry 
for the Ag pseudopotential. In the fully relaxed geometry, 
the Ag atom — in contrast to Li — occupies an asymmet- 
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FIG. 4. Simulated STM images for the honeycomb 
chain-channel model of Li:3xl for (a) filled states and (b) 
empty states. Black circles indicate the projected positions of 
Si and Li atoms in the surface layer. 

ric location with respect to the [TlO] direction and thus 
breaks this mirror-plane symmetry. With respect to the 
bulk crystal, the asymmetry is 0.66 A, in very good agree- 
ment with the diffraction result of 0.58 A g. Moreover, 
this asymmetry distorts the triangular spacing observed 
for Li:3x 1, and leads instead to a pairing of atoms a and 
d from adjacent chains. Physically, this asymmetry is 
driven by the formation of weak Ag-Si bonds, which fa- 
vors two-fold coordination of the Ag atom — in contrast 
to the purely electrostatic interaction of Li and Si, which 
favors three-fold coordination of the Li atom. Further- 
more, this symmetry breaking is accompanied by a dis- 
tortion of the surface-layer Si atoms in which the chain 
hexagons are rotated by ~7° (these distortions were not 
modeled experimentally). This rotation leads to bond 
strains which can be partially relieved by alternating the 
sense of the Ag asymmetry in adjacent channels. Indeed, 
we find that this 6x1 reconstruction, which has paired 
Si atoms but unrotated hexagons, is lower in energy by 5 
meV/A 2 relative to Ag:3xl. These differences have been 
widely observed in comparative STM studies, which in- 
deed find asymmetrically paired maxima and 6x1 peri- 
odicity for M=Ag, but not for Li 0l|]. 

Finally, we suggest a reinterpretation of a recent pho- 
toemission experiment for Mg:3xl |l(|. Since Mg is di- 
valent, a band description of Mg:3xl must be metallic. 
But high-resolution ARPES data shows no peak crossing 
the Fermi level, and on this basis was interpreted as evi- 
dence for a Mott-Hubbard insulator phase |lq| . Our cal- 
culations for Mg:3xl suggest a different interpretation, 
motivated by two important results: (1) The fully relaxed 
HCC model for Mg:3xl is very similar to Li:3xl; in par- 
ticular, the mirror plane normal to [112] is again nearly 
perfect. Thus, by the symmetry argument made above 
for Li:3xl, the matrix element for odd-parity states will 
vanish when the incident light is polarized parallel to the 
mirror plane. (2) The band structure of Mg:3xl is very 
similar to Li:3xl, with one striking exception: the 
state is pulled down from the conduction manifold to sit 
approximately midgap, pinning the Fermi level at ^0.5 
eV above the valence-band edge. The ARPES data shows 



two surface states (labeled SS[ and SS2 in Fig. 3 of Ref. 
Jl6|) whose dispersion is in excellent agreement with S 1 ^ - 
and 5+ in Fig. 0. The odd-parity S% state is not ob- 
served in the data; we propose that 5j~ is not observed 
for the same reason. In other words, Mg:3xl is likely 
a normal band metal, but appears insulating as a result 
of selection rules. Photocmission using unpolarized light 
would, we predict, show a peak at the Fermi level. 

In summary, we have shown theoretically that 
the honeycomb chain-channel reconstruction model for 
A// Si (111)3x1 is energetically far more favorable than 
any other model proposed to date. The stability of this 
model is largely due to the formation of true Si double 
bonds in the surface layer. For Li:3xl and Mg:3xl, the 
existence of a mirror-plane symmetry implies that the 
surface states have definite parity — consistent with the 
findings of angle-resolved photoemission data. Finally, 
this model gives a complete accounting of both filled- 
and empty-state STM data, including the important dif- 
ferences between Li:3xl and Ag:3xl surfaces. 
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